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In this paper we present a theoretical model for the dewetting of ultra-thin polymer films. Assum- 
ing that the shear-thinning properties of these films can be described by a Cross-type constitutive 
equation, we analyze the front morphology of the dewetting film, and characterize the time evolution 
of the dry region radius, and of the rim height. Different regimes of growth are expected, depending 
on the initial film thickness, and on the power-law index involved in the constitutive equation. In 
the thin-films regime, the dry radius and the rim height obey power-law time dependences. We then 
compare our predictions with the experimental results obtained by Debregeas et al. [Phys. Rev. 
Lett. 75, 3886 (1995)] and by Reiter [Phys. Rev. Lett. 87, 186101 (2001)]. 

PACS numbers: 68.60.-p, 68.15.+e, 68.55.-a, 83.10.-y 
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I. INTRODUCTION 

Thin liquid films have considerable scientifical and 
technological importance, and have numerous applica- 
tions. In engineering, for instance, they serve to protect 
surfaces, and applications arise in paints, adhesives and 
membranes [jlj . Thin liquids films display a variety of in- 
teresting dynamics phenomena and have therefore been 
the focus of many experimental and theoretical studies 
0. In particular, thin polymer films have recently at- 
tracted a lot of interest since understanding their prop- 
erties such as viscosity Q , chain mobility Q and stability 
H is essential for optimization. 

A long time ago, Taylor |6| and Culick JtJ analyzed the 
growth of a circular hole in a thin liquid sheet M . By 
balancing surface tension forces against inertia B, they 
found that the rim of liquid at the edge of the films re- 
tracts at a constant velocity, a prediction first checked ex- 
perimentally by Mc Entee and Mysels |h], [ll| . Debregeas 
and collaborators have recently studied the bursting 
of thin suspended films of very viscous liquids. These 
experiments revealed unexpected features: (a) First, the 
retraction velocity grows exponentially with time (with 
a characteristic time scale Tj = hirj/\S\, where hi, r] and 
S are respectively the initial film thickness, the viscos- 
ity and the spreading coefficient (b) second, the 
liquid is not collected into a rim and the film remains 
flat through the retraction. According to these authors, 
the uniform thickening of the retracting film was a con- 
sequence of its viscoelasticity, which permits an elastic 
propagation into the film of the surface tension forces 
acting on the edge. Brenner and Gueyffier Q showed, 
however, that the absence of rim can also result from a 
purely viscous effect jlj| [l(J . 

Very recently Jl7[j , Reiter studied the dewetting of ul- 
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trathin (i.e. thinner than the coil size), almost glassy 
polystyrene (PS) films deposited onto silicon wafers 
coated with a polydimethylsiloxane (PDMS) monolayer. 
He found that a highly asymmetric rim, with an ex- 
tremely steep side towards the interior of the hole and 
a much slower decay on the rear side, builds up pro- 
gressively [|l8] , with the maximum height increasing lin- 
early with the diameter of the hole. We recently E^] 
proposed a theoretical model to explain such deviations 
from the behavior of pure liquids, based on the shear- 
thinning properties of the polymer film. Assuming that 
the stresses inside the film saturate logarithmically with 
the strain rates, we showed that different regimes of 
growth are expected, depending on the initial film thick- 
ness and the experimental time range. Other theoreti- 
cal approaches for polymer film dewetting have recently 
been proposed by Herminghaus et al. pc| and Shenoy et 

In this paper we aim at precising our methods of res- 
olution, and characterizing the dewetting process for a 
Cross-type constitutive equation (interpolation between 
a viscous behavior at low strain-rates and a power-law de- 
pendence for the viscosity versus the shear strain rate), 
following recent experimental results by Dalnoki-Veress 
et al. for PS films |P2| . We study the profile of the dewet- 
ting film, and characterize the time evolution of the dry 
radius, and of the rim height (as summarized in appendix 
^). Finally, we discuss the adequation between our theo- 
retical predictions and experimental results by Debregeas 
and Reiter. 



II. THE MODEL 

Figure [j] shows the film geometry. h(r, t) is the profile 
of the film, h m (t) is the height of the rim, and Rd(t) is the 
radius of the dry zone. vOr, t) is the radial, axisymmetric 
flow field (cf. appendix Eh. On a non-wettable, smooth 
and passive solid substrate like the PDMS-coated silicon 
wafer used by Reiter, this plug-flow description is valid as 
long as hi <C b where b is the hydrodynamic extrapolation 
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the power-law index n strictly between and 1. At low 
strain-rates (jr < A; -1 /™), this law displays a viscous- 
type behavior (cr m ss 7707, with a zero-shear viscosity 
Vo = vqt), while for large values of jr, a m obeys the 
well-known power-law model. This popular expression, 
which will be used all along with this article, is applica- 
ble to a number of polymer materials and complex fluids 



FIG. 1: Film geometry: h(r,t) is the profile of the film, h m (t) 
is the height of the rim, and Rd(t) is the radius of the dry 
zone. The initial (t = 0) step-like profile is represented by the 
dashed line. v(r, t) is the radial, axisymmetric flow field. 



B. Equations 

Assuming the fluid to be incompressible, mass conser- 
vation leads to: 



length (cf. |2q]). The range of thicknesses studied by 
Rcitcr allows this simplifying assumption. 



Constitutive law of the material 



In order to characterize the rheologic properties of the 
film, we introduce the stress tensor cr^- = —pSij + cr™j, 
where p is the pressure and <7"\j represents the effects of 
internal friction. We relate <r m i4 to the strain rate tensor 
jij by a constitutive law of the form: 



cr m (7) = a mr), 



(1) 



where <tq and t are material constants, and $ is a 
generic function. 

For a purely viscous liquid, the function <!> is linear: 



4>(7t) = jr 



(2) 



In this case, the fluid viscosity rj = a" 1 /^ is constant 
and equals the zero-shear viscosity r/o = oqt. 

For polymers, just above T g , it is expected within the 
framework of the free-volume model that a m is expected 
to vary logarithmically with 7 as E4] : 



$(77-) = ln(l + 77-) 



(3) 



The features of the dewetting regimes obtained with 
this rheological law were presented elsewhere [Q, and 
will only be briefly mentioned in the present paper in 
Appendix [§]. As mentioned in jl9| , this logarithmic law 
does not permit a simple analytical description of all the 
regimes of time and thicknesses. 

An alternative way of describing the polymer rheology 
is to use the well-known Cross model |E5|: 



$(7r) 



1 + fc(7r) n ' 



(4) 



1 dh(r, t) r/3(r, t) dh(r,t) 
:( KT~ + Z fr—) 



a(r,t) - p{r,t) 



h(r, t) y dt 



(5) 



Equation ^| involves two unknown functions air, t) 
and f3{r,t), that are positive, dimensionless forms of the 
strain-rate components 7 rr and 7, 



a = — T7 r 



8 T 



(6) 



We thus need two additional equations to determine 
h(r,t). First, note that the following partial differential 
equation can be directly derived from Eq.^: 



dp 
dr 



(3 + a 



(7) 



Neglecting the inertial term, conservation of momen- 
tum (projected on the radial direction (27|) leads to: 



da rr 
dr 



= 



(8) 



It can be shown that Eqs. and ||, along with the free- 
surface boundary condition (i.e. <r zz = at the contact 
with ambient atmosphere pqj), allow one to express the 
strain rate a as a function of (3 only : a = F(J3). Sub- 
stitution of the shear-thinning constitutive law |l] in Eq.|8| 
then leads to the following differential equation for F: 



d[$(FG8)) - $(F(/3) - /?)] _ $(F09)) - $(/3) 



dp 



F(p) + P 



(9) 



k and n being dimensionless constant parameters. A 
shear-thinning behavior is taken into account by taking and P 3> 1, as explained below (cf.IIIAl) 



which has to be solved along with the condition F(0) = 
(far away from the perturbed central region, the strain 
rates must decrease to zero). In FigJ^, we present the 
form of functions a — F(P) for the Cross model, taking 
for each curve k = 10 and different values of power-law in- 
dex n between and 1 . This function can be approached 
by simpler expressions in the two opposite regimes P <C 1 



FIG. 2: The radial strain-rate, a, and the orthoradial strain- 
rate, 3, are related by a = F{0) (see Eq.^). In the purely 
Newtonian case ($(-yr) = yr), the function F is linear: 
a = 3. The other curves represent F for the Cross model 
(EqMj), taking different values of n for a given value of k — 10: 

n 1/4; 1/3; 1/2; 2/3 (an important C3jS6 ; £LS precised in part 

IV B|); and n = 3/4. The two limit regimes /3 < 1 and /3 > 1, 



and the corresponding approached expressions for F are also 
represented. 



In order to solve Eqj^, we should supply it with initial 
and boundary conditions. Our initial profile is assumed 
to be uniform, with a thickness hi, except in a bored 
region ranging from r = to r = Rq (as shown in fig.0). 
Ro is our characteristic radial length used thereafter to 
make r dimensionless: 
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FIG. 3: Variation with the power-law index n of the param- 
eter V( n )- Some other quantities of interest are also repre- 
sented (see text). 



Thereafter, all thicknesses will be made dimensionless 
by normalizing with a characteristic length h* = \S\/oq. 
For Reiter's experimental conditions, at T w 105°C, we 
estimate h* to be of the order of 500 A. This parame- 
ter h* is of crucial importance in our model, because it 
discriminates two regimes of growth whose features are 
quite different. The films whose initial thickness hi is 
larger than h* will be considered as "thick" films, while 
the others (hi <C h*) will be the "thin" ones. 



hi t — 0) — i ^ ^ r — -^o = 1 /iq\ 
^ ' ' 1 otherwise ^ ' 

We do not consider here the origin of the initial dewet- 
ting process: experimentally, it is found that a thick PS 
film on a silicon substrate is metastable and dewets via 
nucleation and growth of dry patches j2{|, while thin- 
ner films (hi <100 nm) are unstable and dewet by spin- 
odal decomposition |30| |3l[] . In his latest experiments on 
ultra-thin films JItJ , Reiter characterized the early stage 
of the dewetting process by the formation and coalescence 
of little holes, with the displaced material uniformly dis- 
tributed between the holes, without visible rims. Our 
initial time t — might correspond to the end of this 
preliminary process. 

Equation || applies outside the dewetted region 
(Rd(t) < r < oo). At the edge of the rim (r = Rd(t)), 
the rim height, h m (t), can be determined by taking into 
account capillary forces. The applied force on the rim, 
pushing the film away the dry area, must be balanced 
by the internal radial stress: \S\ = \a rr \h m (t). Assum- 
ing the lateral extension of the film to be large enough, 
the film thickness must reach hi far from the dry region: 
linv^oo h(r, t) — hi (Vf). The complete resolution of our 
set of equations |[|^-|| can be achieved using a method of 
characteristics (cf . appendix pi) . 



III. RESULTS 

Note first that for a purely viscous liquid ($(-yr) = -yr), 
our model leads to a constant and uniform thickness for 
the film, with an exponential growth of the dry radius: 



Rd(t) = Roe°oX = Roe^i 
h m (t) w hi 



This is in complete agreement with the experimental 
results of Debregeas et al. for the dewetting of suspended 
polymer films (l2| and supported films |3^] . 

How are affected the growth laws of the dry radius 
Rd(t) and the rim height h m (t) by the shear-thinning 
properties of the material ? 

In order to discuss the time dependence of the rim 
height and the dry radius, we need to compare hi with 
the characteristic thickness h* . 
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The rim height obeys the equation: 



0.2 0.4 



T 



FIG. 4: Short-time evolution of the rim height h m : plot of h m 
versus t/r, for hi = 0.1, 0.5, 1 and 2 h* . Note that the curves 
corresponding to hi =0.1 and 0.5h* cannot be distinguished 
within the graph scale. 



A. Thin films regime (hi <C h*) 

1. Initial stages of growth 

The range hi <C h* allows different regimes of growth 
for the initial stage of hole formation, depending on the 
specific characteristics of the constitutive law. 

Simple analytical laws can be obtained with the Cross 
law ($(7t) = (tt)/(1 + /c(7t)™)). It can be shown that 
the solution a[/3] of EqJ| is characterized by a ~ j3 for 
low strain-rates, and a ~ ip(n)P for p 3> 1. It is worth 
noting that the coefficient ^(n) does not depend on k, 
but depends on n as a solution of: 

+ _ l) 1 "™)^ + 1)(1 - n) = ^Z- 1 -" + 1 (12) 

As shown in Fig.||, this coefficient ranges from ip — 1 (for 
n — > 0) to tp — > +oo (for n — > 1), and is involved in the 
laws of growth relative to the rim height and dry radius, 
as shown below. 

The strain-rates a m (t) and P m (t) at the rim edge (i.e. 
for r — Rd{t)) are simply related to the rim height by 
balancing the radial stress acting on the rim and the cap- 
illary forces. This gives the following equation: 

h C^m . Pm 

~Kr ~ 1 + fca" 1 + k(a m - p m ) n ' 



(13) 



admitting the approached forms p m ~ h*/h m for large 
values of h m /h* , and /3 m ~ /i(n, fc)(/i*//i m ) 1 / 1 ~™ for the 
thin films regime. /^t(n, k) is a positive dimensionless pa- 
rameter depending on n and k as: 



fc 



(14) 



V'(n) 1 -" + (^(n) - I) 1 "™ 

These expressions of strain-rates at the edge of the rim 
versus the rim height enable us to find analytical ex- 
pressions of the growth laws during the initial stages of 
dewetting for thin films. 



1 Dh m 



Pn 



(15) 



where D 



di 



stands for the particular derivative. 



In Fig.|4j we present our results about the short-time evo- 
lution of rim height for different initial thicknesses, with 
t ranging from to r. The solution of Eq.[l5] along with 
the initial condition h m (t = 0) = hi is: 



h rn (t) = hi ( 1 + A(n, k, hi)- 



(16) 



with X(n,k,h l ) ee k jl~l 1 (^-) 1 / 1 - n . The adequation 
between our numerical results and this theoretical pre- 
diction is excellent. 

The dry radius is given by: 



Rd 

Rd 



Pm(t) 



and we thus obtain: 



R d (t) =Rq[1 



t \ * 
X(n, k, hi)~ 



(17) 



(18) 



It is physically understandable that Rd, at a given time t, 
is larger for thinner films : for a given applied force \S\ per 
unit length, the thicker the film, the more the material 
to be displaced and the lower the dewetting velocity. 

Equations [l6| - |l8] allow one to define a crossover time 
t c = r/X(n,k,hi). At short times (i.e. t <g; t c ), h m and 
Rd vary linearly with time: h m ps hi(l + (1 — n)t/t c ) and 
R d w i?o(l + %z^V*c)- At longer times (t > t c ), the 
time dependence of h m and Rd is a power-law: h m ~ 
and R d ~ tC 1 -")/^- 1 ). Note that the thinner the film, 
the smaller the crossover time t c . 

A simple expression of Rd versus h r , 
derived from Eqsjl| and || 



can be directly 



Rd(t) — Rq 



h m (t) 



(19) 



This relationship between the dry radius and the rim 
height is shown in Figj^ for different values of power-law 
index n. Note that, for n = 2/3, the expression jl9| ) 
turns out to be linear: this important property w ill be 
discussed in the light of Reiter's experiments in part [V B . 



2. Long-time regime of growth 

After the initial stages of hole growth, the driving force 
of the phenomenon, |S|, is distributed over a rim height 
larger than the characteristic thickness h* . The result- 
ing radial constraint becomes weak, and gives rise to very 
small values of strain-rates a and p. Then, to solve Eq.|l5| 
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FIG. 5: For thin films (hi <C h*), during the initial stages 
of growth, the dry radius is related to the rim height by a 
simple power-law depending on the material rheology (Eq.h9|). 
This figure shows the different regimes depending on n: for 
n = 2/3, the dependence is simply linear; for n > 2/3, the 
radius grows slower than the rim height, while for n < 2/3, 
Rd grows fastly compared to the rim height increasing. 
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FIG. 6: Long-time evolution of the rim height h m , for different 
initial thicknesses hi (0.1; 1; 10; 50; 100ft*), with Cross law 
(n = 2/3, k — 5): plot versus (t/r) 1+n to show the linear 
behavior of h m with (t/r) 1 + n after the transient time to. For 
a thick film (hi S> h*), two regimes can be distinguished. 
For t < to ~ t (t^") 1+ "' hm i s constant and Rd(t) increases 
exponentially with time. For t > to, h m growths linearly with 

(t/r) !+" . Note that the curves corresponding to hi =0.1 and 
\h* cannot be distinguished within the graph scale. 



in this regime, we have to find the smallest-order correc- 
tion to the linear behavior of function F: 



(20) 



Including this development in Eq.^, and neglecting terms 
involving second-order (or higher) e, we obtain a non- 
linear first-order differential equation for e. Solving it 
along with the condition e(/3 = 0) = leads to an ex- 
plicit expression for e(/3). Keeping the smallest order in 
development of e for small (3, we find [[33]: 



a ~ (3+ , 

/3<i V n + 1 



■n+l 



(21) 



From Eq.[15| it is then easy to obtain the following 
expression for h m : 



h m (t) = hi ( 1 + v(n, k, hi)- 



(22) 



with v(n, k, hi) = nk(j-) 1+n . We have checked numeri- 
cally (cf. Fig||) that this power-law behavior is obeyed 
in the long-time regime of growth. 

The corresponding law for dry radius in this regime is: 



1 + n ( h 



Rait) = R exp I ^(1 + u(n, k, fc)-)^ - l 

(23) 

For t > T(hi/h*) 1+n /(nk), this law turns out to be a 
stretched exponential: 



R d (t)~R Q em^r^w^ 



(24) 



From Eqjlg, it is easy to d efine a crossover time t e 
between the initial stag es (part [II A 1 ) and the long-time 
regime of growth (part [II A 2 ): 



hi 



(25) 



3. Film profile 

An example of profile obtained for hi = O.l/i* is shown 
in Fig.f^. This profile is characterized by a highly asym- 
metric shape for the rim. There is a striking similar- 
ity between such a profile and those observed by Reiter 
for ultra-thin films |l7| . On fig.|§| we present the time 
evolution of this profile (for t ranging from t = to 
t = 5.10~ 2 r, with a time step of 10~ 2 r). 

For the sake of conciseness, we focus there on the case 
corresponding to the range of thicknesses and times cov- 
ered by Reiter's experiments, i.e. the initial times of 
growth for thin films (cf. |l7) and f_V). For clarity rea- 
sons, let us write fl(t) = 1 + \(n, k, hi)-. In the regime 
hi <C h* , the partial derivative equation [| for h(r,t) re- 
duces to: 



dh(r,t) k _ . 2- n - n i dh , , _ , k . 2- n - n * 



(26) 

The resolution of Eq,26 can be achieved using a 
method of characteristics, as explained in appendix 
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FIG. 7: Film profile for hi = O.lh* and t = 2.1Cn 2 T. The con- 
stitutive equation of the material is the Cross law (Eq.^|) , with 
a power- law index n = 2/3. For clarity reasons, a sectional 
view is displayed on a quarter of circle, and cartesian co- 
ordinates (x,y,z), naturally associated with cylindrical ones 
(r, $, z), are used. 
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B. Thick films regime (hi 3> h*) 

For thick films, as shown in FigJ^ (for instance in the 
case hi = 100ft-*) two different time regimes can be dis- 
tinguished: at the beginning of hole formation, h m is 
nearly constant and equals its initial value hi during a 
long period of time, before growing faster. Here again, a 
distinction between different time regimes is necessary. 

1. Initial stages of growth 

During the early stage, the strain rates at the rim are 
constant and small, as a consequence of the constant 
large thickness. In this case, the rheological law of the 
film is viscous-type and the dry radius increases exponen- 
tially, as for the case previously discussed (see EqJTl]): 

, . si* J_ 

R d (t) = R e'o^, = R e^ ^ 
h m (t) ~ hi 

This viscous-type behavior remains valid up to a 
crossover time to « T {j^) 1+n ■ 

2. Long-time regime of growth 



0.25 - 




FIG. 8: Evolution of a film (Cross model, n = 2/3) profile 
h(r,ti), for hi — O.lh* and ti ranging from t = (initial 
configuration, in dashed lines) to t = 5.1CP 2 r with a time 
step of 10 -2 r. 



The solution is: 



We have checked numerically (cf. Fig.^) that the an- 
ticipated analytical behavior 




is obeyed after the crossover time to- 

Similarly, the preceding exponential law for the dry 
radius connects with the stretched exponential, whose 
power depends on power-law index n (cf. Eq.|24|) : 

R d (t)^Roe^- in+ ^)^ 

Fig.p which present the time evolution of R d for different 
initial thicknesses hi with the Cross model, shows that 
this anticipated behavior is obeyed for t^>to- 

We emphasize the fact that the crossover from a sim- 
ple exponential to the stretched exponential regime is a 
consequence of the non-linearity of rheological law nl 



h(r, t) = hi — (27) 

(l + r V+i _ jRrf (t)V'+l)lf+T 

We checked (cf. appendix [5]) that this expression has 
a wide range of validity: during all the initial growth, the 
strain-rates remain high (a and /3 3> 1) far away the hole 
periphery, and thus Eq.Ej is an acceptable description of 
the whole film profile. 



IV. DISCUSSION 

A. Debregeas' experiments 

G. Debregeas et al. JlJ] carried out experiments on 
viscous bursting of freely suspended films of long-chain 
polymers. The high molecular weights of polydimethyl- 
siloxane used in these experiments lead to high viscosities 
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FIG. 9: Logarithmic plot of the dry zone radius Rd versus 
[t/ T ) n l (!+"■), for hi = 0.1, 1, 10 and 100 h* , in the case of 
Cross law with n = 2/3 and k = 5. Note the linear behavior of 
log 10 [Rd(t)/Ro] with (i/r) n/(1+n) after a transient time: the 
anticipated stretched exponential (Eqj2^) is indeed obeyed. 
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FIG. 10: Velocity field v(r,ti) for different initial thicknesses 
(hi — 0.1; 1; 10/i*), at a given time All radial distances are 
normalized by the dry radius at time U for each thickness. 
A curve 1/x is represented for comparison.lt is clear that the 
velocity field is almost proportional to 1/r for a thick pro- 
file (hi = 10/i*), similar to the long range radial plug flow 
observed by Debregeas. 



(Vo > 600000 cP) bursting processes where viscous dissi- 
pation dominates inertia. 

The range of thicknesses experimentally studied by 
Debregeas et al. (5 to 250 /xm) covers the domain 
hi 3> h*. As the corresponding crossover time is very 
large (to 3> r), these experiments mainly covered the 
first regime t <C to (the criterion hi 3> h* is in fact a ba- 
sic hypothesis of their "soft balloon" model |}4| ). In our 
model, the velocity field v(r, t) is found to be almost pro- 
portional to 1 jr at any time for thick profiles (cf . Figjl0[) , 
similar to the experimentally observed long range radial 
plug flow @. 

Thus, our model accounts well for the exponential 
growth and the absence of rim characteristics of the 
dewetting regime observed by Debregeas et al. for vis- 



cous polymer films. 



B. Reiter's experiments 

G. Reiter jlTj studied the dewetting of ultra-thin, 
almost glassy polystyrene films deposited onto silicon 
wafers coated with a polydimcthylsiloxanc monolayer. 
The thicknesses of the PS films used range from hi =10 
to hi =60 nm. Compared with the characteristic scale 
h* m 50 nm, the ratio hi/h* ranges from about 0.2 to 
unity. The characteristic time r for PS under the experi- 
mental conditions is very long, up to more than a year for 
the lowest temperatures conditions (T ranges from 103°C 
to 130°C, to be compared with the bulk glass transition 
temperature: T g — 97.5°C). We can thus expect that the 
relevant regime for comparing our results with Reiter's 
exp eriment al findings is the limit hi <Si h* , and t/r small 
(cf. |EII Alj) . 

Three years ago, Dalnoki-Veress et al. made an ex- 
perimental study of hole formation and growth in freely 
standing PS films p2fl . In agreement with Debregeas' 
experiments, they observed exponential growth of the 
hole radius and uniform thickening, but with a deviance 
from this regime for long-time regimes of growth, a fea- 
ture they analyzed in terms of shear-thinning properties. 
They indeed characterized the polymer rheology at the 
temperature of their experiments (T = 115°C) and found 
that the best fit for the viscosity 77 versus the shear strain 
rate 7 is a power-law dependence: 



V ~ 171 



(30) 



with a power-law index d = 0.65 ± 0.03. For the fol- 
lowing, let us suppose that the rheology of ultra-thin PS 
films used by Reiter can be described by the Cross model 
(Eq.^j) with such a power-law index: n = 2/3. It is of 
course an hypothesis that would need an experimental 
validation. In particular, the rheology of PS films at 
temperatures very close to the glass transition may be 
quite different. Keeping this point in mind, the follow- 
ing remarks are tentative explanations of some regimes 
of dewetting experimentally observed by Reiter. 

If n = 2/3, Eq.[l2| leads to the exact value ip = 2. Then 
Eq.rLS turns out to be simply linear: 



Rdjt) = h m (t) 

Rn hi 



(31) 



This property was indeed observed by Reiter: the max- 
imum height of the rim grows at the same speed as the ra- 
dius of the hole |l7) . The special role played by this value 
n = 2/2> was also pointed out by Shenoy and Sharma 
in their study of the dewetting of power hardening vis- 
coplastic solid pl| . 

We can also focus on the individual dependences of 
h m and Rd versus time. As precised in part III A 1 , for 
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t 3> t c , the rim height and the hole radius follow a power- 
law time dependence: h m ~ i 1 "™ and R d ~ i(i-»)/M-i). 
In particular, for n — 2/3, we obtain: 



Rd 



(</r)V3 



(32) 



Very recently, Reiter [M characterized the time de- 
pendences of Rd and h m for very thin films. He ob- 
tained several dewetting regimes, and work is in progress 
to analyze these results. It seems that several experimen- 
tal factors must be taken into account for this analysis. 
Among them, the influence of film preparation appears to 
be decisive: the spin-cast PS film is a highly metastable 
form of matter j^]. The sample preparation (including 
the time of healing) can lead to thin films where resid- 
ual constraints on the chains persist, and may have an 
influence on the dewetting process. 

Quite interestingly, despite all these remarks, for a 20 
nm PS film deposited on a 10 nm PDMS layer around 
T = 120°C, the results of Reiter's experiments seem to 
indicate that the maximum height of rim and the hole 
radius obey power-laws: 



Rd 



(t/r) b 



with a ' 
with b ; 



: 0.38 ±0.02 
0.34 ±0.04 



(33) 



The time range on which Reiter calculated these 
power-law indexes is up to about 100 minutes. As for 
a 20 nm film, the ratio hi/h* is smaller than unity, the 
crossover time t c is very short, and we can expect that 
the power-law regime (Eq.|32]) is established. Thus we can 
compare Eqs.[32] and |33| and remark that the two experi- 
mental coefficients are indeed very close to the predicted 
one: 1/3. This adequation between our model and Re- 
iter's experimental results appears to be satisfying. But 
many questions remain open, as the role played by the 
thin PDMS monolayer on which the PS film dewets is 
only taken into account by the free slippage condition 
p6| . In fact, it was shown by several studies that the 
lower layer is deformed by the moving upper one | p7| , 
which can leads to viscous dissipation and to a marked 
decrease in the dewetting velocity Q . 

Let us conclude by a tentative idea: observing by Tap- 
pingMode atomic force microscopy (TM-AFM) the dry 
zone on the PDMS substrate after the removal of PS film, 
Reiter j3i| discovered the formation of straight lines, sev- 
eral microns long, starting to appear at hole diameters 
larger than about 1 /zm. These lines are PS microfibrils, 
whose width can be as low as about 2 nm (to be com- 
pared with the radius of gyration of the polymer in the 
bulk, of order 50 nm). This observation suggests that the 
radial strain-rates are much larger than the orthoradial 
ones at the hole periphery and involve this elongational 
deformation of the molecules. This property could be re- 
lated to the fact that, in our model, the radial strain-rate 
(a) is large compared with the orthoradial one (/3) dur- 
ing the initial stages of growth for thin films (cf. FigJ3, 



regime |3 > 1). In particular, note that this dissymme- 
try between values of a and j3 is a direct consequence of 
the shear-thinning behavior of the polymer: for a purely 
viscous fluid, a — (3 for any thickness. 



V. CONCLUSION 

In conclusion, our model accounts well for the expo- 
nential growth and the absence of rim characteristics of 
the dewetting regime observed by Debregeas et al. for 
viscous polymer films. Taking into account the shear- 
thinning behavior of the polymers near T g enables us to 
see the modifications induced by this particular rheol- 
ogy on the film morphology. It appears from our results 
that the early stages of dewetting for a shear-thinning 
polymer film are mainly determined by the ratio of its 
initial thickness hi to a characteristic scale h* (related to 
the driving force of the process, S, and the rheological 
response of the material, characterized by <7o). 

The film profile exhibits a sharp, asymmetric rim simi- 
lar to the one observed by Reiter in his latest experiments 
with ultrathin PS films. The rich variety of growth be- 
havior is summarized in appendix In all cases, the 
long-time evolution of h m is a power-law (whose coeffi- 
cient depends on the rheology of the material studied), 
and Rd obeys a stretched exponential law. When taking 
for the power-law index n the value obtained by Dalnoki- 
Veress et al., in the thin- film limit, our predictions for the 
dry radius and the rim height time-evolutions seem to ex- 
plain some of Reiter's latest experiments. In particular, 
the linear dependence of the rim height versus the dry ra- 
dius experimentally observed is predicted by our model 
for the power-law index n — 2/3. 

Work is now in progress to incorporate in our model the 
Laplace pressure (see pal ) which might lead to oscillatory 
dewetting fronts 14, ]lq|, and to evaluate the possible 
effects of chain entanglements. In the future, we aim to 



study the dewetting of thin polymer films below T g , where 
the existence of a yield stress in the rheological response 
of the material may lead to new dewetting morphologies. 



Acknowledgments 

We are grateful to G. Reiter for very useful discus- 
sions and for providing us with his experimental results 
on the dewetting of ultra-thin polystyrene films prior to 
publication. 

We also wish to acknowledge fruitful discussions with 
A. Aradian, G. Debregeas, and M. Sferrazza. 



APPENDIX A: A FEW REMARKS ABOUT THE 
VELOCITY FIELD INSIDE THE FILM 

Our choice of model (see |l|) for the dewetting film 
deserves a few comments. As defined above, v(r, t) is 
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the radial part of flow field. The incompressibility of the 
material implies the additional existence of a small, but 
non-zero vertical component of the velocity field: let us 
note w this vertical part. 

We already argued that the characteristic thickness of 
films studied by Reiter and Debregeas allows a plug-flow 
assumption for the film: it implies that the radial part 
v(r, t) is independent of variable z. This hypothesis does 
not imply the same property for w: we must keep for it 
the most general dependence w(r, z, t). 

In cylindrical coordinates, for a velocity field of the 
form (v r ,v$,v z ), the liquid incompressibility leads to: 



dv r 
dr 



V 



1 dv$ 
r <9$ 



dv z 
dz 



= 



Here, as i>$ = for symmetry reasons, one has: 



dv 

dr 



dw 
dz 



(Al) 



(A2) 



Ind epe ndence of v versus z allows a simple integration 
of Eq. A2 at a given r: 



lz=0 

which gives: 



(fl- + -)<k 

or r 



h{r,t) 



z=0 



dw 
dz 



dz, 



(A3) 



h(r, t).(^- + -) = -w (r, z = h(r, t),t) 
dr r 



(A4) 



The vertical velocity at the top of the film, w(r, z — 
h(r,t),t), can be easily related to the particular deriva- 
tive of h(r,t), provided that the film thickness remains 
small compared with the characteristic lengths of vari- 
ation of the radial speed component and film thickness 



/ , / % \ ,dh dh, ... 

w (r,z = h(r,t),t)n(— + v.—) (A5) 

Thus, in a more rigorous way, we recover the incom- 
pressibility condition (see Eq.||): 

1 dh{r,t) rP(r,t)dh(r,t) a(r,t)-/3(r,t) 



h(r, t) v dt 



dr 



The plug-flow assum ptio n has a direct consequence: 
the left-hand side of Eq.jXJ is independent of z; so does 
the right-hand side, i.e. dw/dz, which implies that w 
has a simple linear dependence versus z: 



w(r, z, t) — w(r, z = h(r, t), t) 



h(r,t) 



(A6) 



To finish with this point, note that taking into account 
this small but non-zero vertical component of the velocity 



field induces the presence of a supplementary term in the 
conservation of momentum (equation ||) : 



da r 



dr 



do~ rz <7 r 
dz 



Q 



(A7) 



The possible consequences of this supplementary fric- 
tion force on our predictions for dewetting regimes will 
be studied in a forthcoming publication. 



APPENDIX B: 



RESULTS FOR LOGARITHMIC 
FORM OF $ 



For polymers, just above T g , it is expected within the 
framework of the free-volume model that a m varies log- 
arithmically with 7 as K4] : 



<&(7t) = ln(l + jt) 



(Bl) 



At low strain-rates (tt < 1), this law displays a 
viscous-type behavior (o~ m ~ ?7o7: with a zero-shear vis- 
cosity t]q = o~qt), while for large values of jr, o~ m reaches 
an almost constant value (shear-thinning behavior). 

In order to deal with a minimum number of parame- 
ters, we assume that o~ m (j) is an odd function (c m (7) = 
-cr ln(l + |t|t) if 7 < 0). 

This appendix briefly presents the main results for thin 
and thick films. 

For thick films, at short times, the dry radius and rim 
height evolve as for a purely viscous fluid (cf. Eq.pT|): 



J R d (t) = Roe^H = J R e T * 
\ h m (t) w hi 

For thick films, at long times (t > to ~ T(hi/h*) 2 ), the 
rim height grows proportionally to \ftfr: 



while the dry radius expands like exp (4-y/i/r): 



Rd{t) ~ R 
i>to 



(B2) 



t/r 



(B3) 



This law is also a stretched exponential, as for the 
Cross law. This property is due to the fact that the 
function a — F{j3) is analytical (in the mathematical 
sense), i.e. admits a series expansion near whose first 
non-linear term is 1 . It indeed corresponds the limit 
behavior obtained by taking the Cross law with n — ► 1. 

For thin films, at short times, a sharp rim grows very 
quickly (with an initial speed h m ~ hi exp h*/2hi). In 
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this regime, h m (t) is well described by the integral equa- 
tion: 



dx ~ — 

x *<t T 



(B4) 



Along a given characteristic, the rim height h^(t) obeys 
a simple equation: 

dt r e (i)^ +1 tw 
whose solution is explicit and of the form: 



For thin films, at long times, as for the Cross law, the 
hole growt h conn ects with the long-time regime of thick 
films (Eqs.pgpf). 



APPENDIX C: SUMMARY OF THE 
DEWETTING REGIMES FOR A CROSS 
CONSTITUTIVE EQUATION 

Figures [ll] and [l2] propose a summary of the differ- 
ent dewetting time regimes for a Cross constitutive law, 
depending on the initial thickness (to be compared with 
h*), and relative to different time ranges (to be compared 
with crossover times t c , t e and to). 



Inverting Eq.D2 to express £ as a function of r and t 
enables us to extract the complete equation of the profile 
(Eqf7l): 



h(r, t) = hi 



Range of validity 



APPENDIX D: ANALYSIS OF THE FILM 
PROFILE 

Method of characteristics 

To know the profile h(r, t), the equation ^6] has to be 
solved: 



dh(r,t) k 
dt + 



n(ty 



i* dh 
Or 



al ong with the precedently precised initial condition 
(Eq.|lO) and the boundary values h(r — > oo, t) — hi (Vt) 
and h(r = Rd{t),t) = fe m (t) . As a a quasi-linear partial 
derivative equation, Eq.g6| can be solved using a method 
of characteristics, as precised in this appendix. 

Writing dh(r, t) — ^rdt + dh g^ t ' > dr, we can replace the 
partial derivative with respect to t thanks to EqJ|^. Let 
define a characteristic by the function r^(t), indexed by 
the parameter £ with rg(i = 0) = £, and obeying the 
following ordinary differential equation: 



dt 



n(ty 



(Dl) 



Eq.Dl is tractable, and gives each characteristic as a 
function of time: 



■<i)+i 



n(ty 



i 



(D2) 



In Fig.[T| is shown a representation of characteristics 
r(:(t) versus time, for a thin film (hi — 0.1/i*). Note that 
these curves were extracted numerically from our model, 



and are not a graphic illustration of solutions of Eq. D2 



In the film zone near the hole (r > Rd(t)), the strain- 
rates a and (3 are high. Thus Eq.[7] gives: 



P(r,t) =(3 m (t) 



( R d (t) 



n{ty 



(D5) 



Our description of film profile remains valid as long as 
P( r ,t) > 1- Let us define the radial distance r*(t) up 
to which this condition is verified : (3(r*(t),t) — 1. The 
ratio r*(t)/Rd(t) is given by: 



"(*) 



R d {t) 



= nVL(t) 



(D6) 



We regard the equation of profile as valid if r*(t) ^S> 
Rd(t), and arbitrarily choose a criterion r*(t) > TRd(t) 
with a factor T ~ 10 for instance. It is then found that 
r*(t) > TRd(t) as long as time t is smaller than a critical 
value t* given by: 



t* = 



1 - n fi(n, k) ( h* 



(D7) 



For hi = 0.1/i*, n = 2/3, k = 5, and T = 10, Eq.|D7 
shows that the high strain-rates hypothesis remains valid 
on a wide region as long as t < 1000 r. Thus, the Eq.^ is 
a good description of the film profile during all the initial 
regime of growth. 



APPENDIX E: INFLUENCE OF THE INITIAL 
SHAPE OF FILM PROFILE 

The results presented above are obtained taking a dis- 
continuous shape for the initial profile of the film: 

h(r t — 0) — { ^ ^ r — -^o = 1 
^ ' ' I otherwise ^ ' 
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h m (t) = h i (l+^) 1 - n 






1 


(= 



Exact 
solutions 

Crossover 
times 



Approximated h ^ „ t J/j m ( t ) ^ t l-n 
time 



|i? rf (f) = Ro e 



nst.((l+v ilHS-l) 




dependences Rd(t) ~ i M.Rd(i) ~ i i~^"^H-RdW ~ ex P 



const, t 1 ^ 



Case 
n = 2/3 



M*) ~ * |M*) ~ i 1/3 I 



M*) ~ * 3/5 I 



~ * Ii? d (t) ~ t V^M-fidW ~ exp 



ist. i 2/5 



FIG. 11: Summary of the different time regimes for the dewetting of thin films (hi <C ft*). Depending on each time range, 
the exact solutions and correspo nding approximated time-dependences are precised. We also point out the laws for the special 
value n = 2/3 (discussed in part IV B) 



Exact 
solutions 

Crossover 
time 

Approximated 

time 
dependences 




_h m (t) = h 
R d (t) = R e 





h m (t) ~ const 
R d (t) ~ exp [t] 



h m {t) ~ t — I 



mRd(t) ~ exp const. 



FIG. 12: Summary of the different time regimes for the dewetting of thick films (hi ^> ft*). Depending on each time range, the 
exact solutions and corresponding approximated time-dependences are precised. 



Does this discontinuity at r = i?o play a role in the hole 
growth? In other words, is it the source of the marked 
asymmetry observed for the profiles obtained with our 
model? 

To probe this hypothesis, we made the same resolution 
of our system of equations, starting with a continuous 
profile: 



h(r,t = 0) 







if r > Ro 
if r < Ro 



1 



(E2) 



This profile is infinitely derivable at r = Ro, the pa- 
rameter x g ivm g & n order of the characteristic length of 
transition between the regions h = and h — hi- The 



profiles obtained with this initial condition are shown in 
figure |lj, with x = 0-1 (the shape of the initial profile, 
represented in dashed lines in Fig.|l4|, is very smooth). 
The rim growth seems to be very fast for a thin film (the 
time range of curves in Figgis to 2.1CT 2 t). Similarly, 
the increasingly asymmetric shape of the profile shows 
that, even when starting with a smooth profile, the re- 
sulting morphology of the rim is the same as in Rcitcr's 
experiments. These considerations prove that the discon- 
tinuity of the initial profile is not a necessary condition for 
the appearance, growth and stiffening of an asymmetric 
rim as those observed by Reiter. 
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